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The  research  carried  out  with  the  support  of  this 
contract  during  the  first  year  of  its  existence  has  been  con¬ 
cerned  with  the  following  problems: 

(a)  re -determination  of  the  libration  constants 
of  the  Moon; 

(b)  evaluation  and  measurements  of  the  old  Paris 
photographic  plate  of  the  Moon  -  covering  the 
full  range  of  libration  -  for  determining  the 
departures  of  the  actual  surface  of  the  Moon 
from  a  sphere;  and 

(c)  development  of  methods  for  the  determination 
of  small  differences  in  height  (or  depth)  of 
very  shallow  surface  features  -  such  as  the 
lunar  domes,  rills,  and  wrinkle  ridges  - 
above  (or  below)  the  surrounding  landscape. 

Distinct  contributions  to  all  three  of  these  lines  of 
research  have  been  made  in  the  course  of  the  year;  and  the 
progress  made  along  each  will  be  discussed  below  in  turn. 

The  work  under  (a)  was  performed  by  Professor  K.  Kosiel 
of  the  University  of  Krakow,  President  of  the  I.A.U.  Commission 
No.  IT  (SelenodeayX  who  visited  Manchester  for  a  period  of  six 
months  (October  1961  through  Marqh  1962)  largely  to  use  the 
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University's  electronic  computer  for  this  purpose.  In  this 
work  Professor  Koziel  was  joined  by  his  Krakow  associates, 

J.  Mietelski  (October -December  1961)  and  J.  Maslowski 
(January -March  1962),  ably  assisted  at  the  electronic  computer 
by  Mrs.  Mary  Gorman  of  our  own  staff. 

Preparations  for  work  on  task  (b)  have  been  made  in 
collaboration  with  Dr .  Th.  Weimer,  of  the  Paris  Observatory, 
to  modernize  first  the  equipment  available  to  him  at  Paris  for 
the  purpose  of  the  contemplated  research;  while  work  under  (c) 
was  pursued  throughout  the  year  at  Manchester  by  Messrs 
D.  Dale  and  T.  W.  Rackham,  under  the  supervision  of  Professor 
Z.  Kopal.  The  contributions  of  Mr.  T.  W.  Rackham  to  three- 
dimensional  topography  of  the  Moon  earned  him  the  Master's 
degree  at  the  University  of  Manchester  in  December  1961;  and 
it  is  expected  that  Mr.  Dale’ s  efforts  in  this  field  will  earn  him 
the  same  academic  reward  during  the  coming  session. 

(a)  Re-determination  of  the  Libration  Constants 
of  the  Moon. 

The  importance  of  determining  the  exact  shape  of  the 
Moon,  and  of  the  parameters  (such  as  the  difference  between  the 
principal  moments  of  inertia)  characterizing  the  internal  structure 
of  the  lunar  globe  can  hardly  be  overestimated  at  the  present  time, 
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when  plans  are  under  way  to  land  men  on  the  Moon  before  the 
end  of  the  present  decade.  But  the  first  astronauts  landing  on 
the  surface  of  our  satellite  will  need  good  charts  of  the  Moon, 
similar  to  the  charts  of  the  Earth;  and  the  photographs  -  even 
the  best  ones  -  cannot  replace  such  charts.  The  basic  data  for 
the  construction  of  these  charts  must  come  from  selenodetic 
work  -  similar  in  method  to  terrestrial  geodesy,  except  for  the 
fact  that  all  measurements  must  be  carried  out,  by  means  of 
suitable  telescopes,  at  a  distance  which  never  becomes  less 
than  356,000  km.  The  tasks  of  selenodesy  are  further  compli¬ 
cated  by  the  fact  that,  on  the  Moon,  we  do  not  possess  any 
obvious  surface  of  reference  corresponding  to  the  terrestrial 
sea-level.  In  fact,  the  only  way  to  establish  the  deviations  of 
the  actual  surface  of  the  Moon  from  a  sphere  is  from  the  stereo¬ 
scopic  effects  exhibited  by  such  deformations  in  the  course  of 
lunar  librations . 

The  observations  of  apparent  motions  of  specific 
points  of  the  lunar  surface  with  respect  to  the  limb  of  the  Moon 
in  the  course  of  a  full  libration  cycle  have  been  one  of  the 
principal  tasks  of  selenodesy  in  the  past  century;  and,  as  is  well 
known,  the  most  important  instrument  designed  specifically  for 
this  purpose  has  been  Bessel’s  heliometer.  In  heliometric 
measurements  -  which  constitute  so  far  the  vast  majority  of 
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the  libration  measurements  of  the  Moon  -  it  has  been  customary 
to  refer  the  position  of  the  crater  Meeting  A,  situated  near  the 
centre  of  the  Moon’s  apparent  disk  (Fig.l.)  to  the  illuminated 
limb  of  the  Moon.  This  is  done  by  measuring  the  distance  bq 
(s  observatum)  of  this  crater  from  the  limb  in  specified  position 
angles  p,  in  order  to  determine  the  position  of  Mdsting  A  as 
referred  to  the  centre  of  mass  of  the  lunar  ellipsoid. 

From  Bessel’s  time  until  recently,  the  adjustment 
of  the  observations  of  the  Moon’s  libration  was 
done  in  two  stages.  First,  the  auxiliary  unknowns 
of  the  problem,  (i.e.,  the  corrections  to  the 
rectangular  plane  co-ordinates  of  the  crater 

tl 


Mosting  A)  were  found,  and  only  the  second  adjustment  gave  the 

proper  unknowns  of  the  problem  -  i.e.,  the  corrections  to  the 

selenographic  co-ordinates  X0  ;  as  well  as  to  the  Moon’s 

radius  to  the  crater  Mdsting  A  (so-called  h)  and  also  the 

corrections  to  the  mean  inclination  of  the  lunar  equator  to  the 

ecliptic  3  and  to  the  mechanical  ellipticity  of  the  Moon 
r  B(C  -  B) 

j  ~  AtC  -  A)  ’  w^ere  A,  B  and  C  denote  the  principal 
moments  of  inertia.  From  the  very  beginning  this  manner  of 
treatment  created  great  difficulties  in  the  choice  of  weights  of 
the  right-hand  sides  of  the  observation  equations  in  the  second 
adjustment.  These  right-hand  sides  are  not  independent,  as 
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Figure  1. 

Central  part  of  the  apparent  disk  of  the  Moon, 
including  the  craters  Ptolemy  (centre), 
Alphonsus,  (top)  Flammarion  (bottom)  and 
Mdsting  A  (indicated  by  arrow). 
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they  were  obtained  from  the  same  equations  in  the  first  stage. 
But  the  adjustment  of  such  equations  in  the  second  stage  is 
incorrect  from  the  point  of  view  of  the  least- squares  method, 
and  no  choice  of  weights  can  secure  here  a  correct  solution. 
The  reason  for  such  a  state  of  affairs  was  the  difficulty  in 
the  execution  of  the  very  extended  calculations. 

In  1948  Professor  Koziel  published  a  paper  about 
the  Moon's  libration  in  which  he  reported  some  progress 
with  this  problem:  but  in  195b  he  published  for  the  first  time 
the  results  for  a  special  case;  and  in  1961  he  worked  out  in 
Manchester  a  general  solution  of  the  problem,  correct  from 
the  point  of  view  of  the  least-squares  method.  The  correct 
mathematical  form  of  the  final  observation  equations  becomes: 
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The  differential  co-efficients  A,  A  ,  B,  . G,  G  occurring 

in  these  equations  can  be  found  from  the  formulae  obtained  on 
the  basis  of  differential  formulae  of  spherical  poly gonome try . 

The  adjustment  of  the  libration  observations  whose  left-hand 
sides  contain  the  fundamental  unknowns  of  the  problem  and, 
on  the  right-hand  sides,  the  observed  residuals  (bq  -  sc),' 
permits  the  adjustment  to  be  carried  out  in  accordance  with 
the  principles  of  the  least-squares  method.  Moreover,  the 
present  form  of  the  observation  equations  enables  use  to  be 
made  of  even  the  so-called  quite  incomplete  observation  evenings 
i.e.,  such  evenings  as,  in  the  extreme  case,  consist  of  only  one 
measurement  of  the  illuminated  limb  of  the  Moon’s  disc.  A 
further  advantage  of  observation  equations  of  this  type  is  that 
they  allow  the  introduction,  into  their  right-hand  sides,  of  the 
corrections  to  sq  for  the  irregularities  of  the  Moon’s  limb  or 
any  other  corrections  and  then  a  parallel  adjustment  may  be 
carried  out  without  any  change  to  the  left-hand  sides. 

It  is  interesting  that  the  simplest  and  most  rigorous 
mathematical  form  of  the  observation  equations  as  given  above, 
is  also  the  simplest  for  the  electronic  computer  "Mercury". 

A  total  of  3328  observations  taken  during  340  evenings  in  the 
period  1877-1915  have  been  used.  Of  these,  46  observations 
were  rejected,  as  the  sq  -  sc’s  were  four  times  larger  than  the 
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mean  error  of  one  observation  equation.  However,  in  no  case 
did  a  whole  evening  have  to  be  rejected;  since,  if  we  use  the 
type  of  observation  equation  under  consideration,  it  is  possible 
to  utilise  every  observation.  When  Professor  Koziel  came  to 
Manchester  he  had  already  prepared  the  differences  sq  -  sc 
for  the  following  series  of  observations: 


Strasbourg  Series  (1877-1879)  by  K.  Koziel; 

Dorpat  Series  (1884-1885)  by  K.  Koziel; 

Bamberg  Series  (1890-1912)  by  J .  Maslowski;  and 
Kazan  Series  (1910-1915)  by  J.  Mietelski. 

The  sq  -  sc  were  calculated  by  using  the  formulae: 
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which  represent  the  solution  of  a  spherical  pentagon  on  the 
selenographic  sphere. 


Then, 
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where  h'*  is  the  topocentric  radius  for  Mtisting  A.  After 
computing  ir  "  from  the  formula 
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Turning  now  to  the  adjustments  and  starting  from 
Hayn’ s  value  of  the  libration  constants,  Professor  Koziel  has 
divided  the  adjustment  calculations  into  seven  parts. 

1.  Professor  Koziel  gave  to  Mrs.  Gorman,  for  pro¬ 
gramming  for  the" Mercury"  computer,  the  following  formulae 
which  have  recently  been  published  in  English  for  the  first  time 
in  Physics  and  Astronomy  of  the  Moon,  ed.  Z.  Kopal  (Academic 
Press,  London  1962): 
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(Yakowkin  hypothesis:  dR  =  x  +  PMy) 


2.  Since  an  accurate  value  of  if  cannot  be  obtained 
from  the  astronomical  almanacs,  it  is  necessary  to  make  suitable 
corrections.  Professor  Koziel  derived  at  Manchester  the 
following  formula  for  correcting  the  right-hand  sides  of  the 
observation  equations  for  h  ^  : 


rr 


A  h'  lcc.  r"  f  s  i  u  ([5  -  ir  " ) 
A  \  un  fp-rr") 


Atv  An  1  |  "  ^ Ww  X 

A  (3  At.3  I  I  Cofi  Co  X 


There  ^re  three  right-hand  sides:  without  limb  corrections,  with 
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Hayn  and  with  Weimer  charts.  On  November  14th,  1961,  this 
new  formulae  was  given  to  the  programmer  and  new  normal 
equations  have  been  obtained  starting  from  f  =  0.73.  It  has 
also  been  checked  that  the  mean  value  of  the  corrections 
for  the  period  in  question  alters  the  d  \  in  the  right  direction 
and  does  not  change  other  unknowns . 

3.  On  October  18th,  1961,  Mrs.  Gorman  received 
for  programming  the  formulae  obtained  by  Professor  Kosiel 
in  Manchester  for  starting  on  the  other  side  of  the  critical 
value  of  f  =  0.662’. 

1)  new  formulae  for  d\  and  df)  ;  and 

2)  the  changes  of  right-hand  sides  for  fQ  -  0.60: 
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4.  Next,  the  coefficients  of  the  free  libration  in 
longitude  should  be  computed: 
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The  following  new  formulae  for  the  adjustment  with  free  libration 
were  given  to  the  programmer  on  November  28th,  1961: 
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From  the  Kazan  series  first  approximations  of  the  amplitude 
of  free  libration  are: 

A  =  19"  (without  limb  corrections);  and 
A  =  8"  (with  Hayn  charts). 

These  are  preliminary  results. 

The  last  part  of  the  work  concerns  the  formation  of 
the  coefficients  of  the  normal  equations  for  the  solution  of  the 
problem  in  the  second  and  last  approximations.  The  respective 
formulae  were  given  by  K.  Koziel  to  the  programmer  in  December 
1961  and  February  1962.  Instead  of  the  unknown  df,  two  unknowns, 
dfj  and  df^  have  been  used  in  order  to  avoid  the  nonlinearity  of 
the  problem,  caused  by  the  effects  of  a  small  divisor  (i.e.  the 
co-efficient  to  the  term  of  argument  2u  in  the  forced  libration 
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in  longitude).  It  is  interesting  to  note  that  the  capacity  of 
"Mercury"  was  found  too  small  for  carrying  out  such  extended 
calculations  as  the  adjustment  of  the  Bamberg  series  and 
these  had  therefore  to  be  divided  into  sections.  In  each  part 
we  have  six  sets  of  systems  of  normal  equations  with  6-10 
unknowns . 

Given  below  are  the  preliminary  mean  errors  of 
one  observation: 


Bamberg. 

Kazan. 

Objective  diameter 

184  mm 

106  mm 

Focal  length 

2.7  m. 

1.7  m. 

Without  limb  correction 

+  O'. '80 

+  0\'87 

With  Hayn  chart 

+  0'.'45 

+  0'J48 

With  Weimer  chart 

+  O'. '43 

+  0  '*  52 

Professor  Koziel  has  also  prepared  the  programming 
for  obtaining  the  sc  (divided  into  five  parts)  and  a  big  job  has 
been  done  in  checking  this  programme.  If  possible,  he  would 
like  to  continue  the  libration  calculations  and  extend  them  to 
cover  the  period  1841  to  1945  in  order  to  obtain  not  only  the 
classical  libration  constants,  but  also  the  free  libration  constants 
in  inclination  B,  b  and  in  node  C,  c. 

The  actual  solutions  of  the  normal  equations  for  the 
definitive  (second-approximation)  values  of  the  libration 
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constants  were  not  yet  carried  out  before  Professor  Kcsiel  and 
his  colleagues  left  Manchester  for  Krakow  by  the  end  of  the 
period  covered  by  this  report.  As  is  evident  from  the  fore¬ 
going  account,  the  bulk  of  the  time  available  to  the  Polish 
investigators  at  Manchester  has  been  spent  in  the  formulation 
of  the  coefficients  of  the  equations  of  condition  based  on  lunar 
heliometric  observations  from  the  years  1877  -  1915,  and  of 
the  coefficients  of  the  normal  equations  based  upon  them.  This 
called  for  an  amount  of  computation  which  would  have  been 
wholly  impracticable  without  access  to  a  modern  electronic 
computer  -  such  as  the  Ferranti  "Mercury"  of  the  University 
of  Manchester  -  which  are  not  available  in  Poland.  On  the 
other  hand,  the  solution  of  the  normal  equations  themselves 
does  not  call  for  the  use  of  an  electronic  computer,  and  can 
be  effectively  performed  with  the  aid  of  desk-type  computers 
as  are  available  in  Krakow.  The  final  results  of  such  compu¬ 
tations  are  awaited  with  interest,  and  should  be  published  in 
the  near  future. 

(b)  Evaluation  and  Measurement  of  the 
Paris  Lunar  Negatives. 

In  addition  to  a  modern  re-discussion  of  lunar  helio¬ 
metric  work,  performed  by  the  visual  observers  throughout 


the  past  century,  photographic  material  of  potentially  equal 
value  for  the  determination  of  the  shape  of  the  Moon  has  been 
known  to  exist  at  the  Paris  Observatory,  where  several 
thousand  lunar  photographs  taken  with  the  Observatory's  "grand 
equatoreal  coude"  of  60  cm  aperture  and  18  m  focal  length 
were  accumulated  between  1890  and  1910  by  MM  Loewy  and 
Puiseux.  The  main  value  of  this  material  rests  on  the  fact 
that  it  covers  in  time  the  entire  amplitude  of  lunar  libration. 

On  the  other  hand,  although  many  photographs  of  this  series 
are  of  excellent  quality,  their  emulsions  cannot  compete  with 
more  modern  photographic  plates;  and,  above  all,  on  account 
of  their  age,  these  plates  cannot  be  transported  elsewhere  in 
bulk. 

In  order  to  evaluate  their  potential  usefulness  for 
the  present-day  selenodetic  research,  and  to  undertake  such 
measurements  of  them  as  may  be  desirable,  sub -contractual 
arrangements  have  been  made  with  Dr.  Th.  Weimer,  of  the 
Paris  Observatory  (Carte-du-ciel  division)  staff,  vice-president 
of  the  I.A.U.  Commission  No.  17  (Selenodesy)  and  a  well-known 
expert  in  this  field.  Pursuant  to  these  arrangements,  Dr.  Weimer 
visited  the  United  States  in  the  summer  of  1961  for  conferences 
with  Mr.  Hunt  at  GRD,  AFCRC,  and  Dr.  Markowitz  at  the  U.S. 
Naval  Observatory  in  Washington,  concerning  the  instrumental 
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equipment  which  he  would  need  at  Paris  to  embark  on  such 
measurements  of  the  selected  lunar  negatives  as  may  be 
desirable . 

Agreement  was  reached,  and  plans  drafted  for  the 
modernization  of  rooms  and  improvement  of  the  measuring 
engine  on  loan  at  Paris  from  the  U.S.  Naval  Observatory, 
to  be  done  partly  at  Paris  and  partly  in  the  United  States. 

The  terms  of  a  sub-contract  to  finance  these  developments 
from  the  funds  available  to  us  have  been  worked  out,  and 
signed  on  March  7th  1961,  by  Professor  Danjon,  Director 
of  the  Paris  Observatory,  and  by  Professor  Kopal  on  behalf 
of  the  University  of  Manchester.  In  the  meantime.  Dr.  Weimer 
embarked  on  a  selection  of  the  Paris  negatives,  taken  at 
different  angles  of  libration,  on  which  craters  with  dimensions 
between  2-5  kms  are  measurable.  A  list  of  such  objects 
suitable  for  measurement  was  published  by  Dr.  Weimer  in 
an  Appendix  to  the  Proceedings  of  the  First  Conference  on 
Lunar  Topography,  held  at  Bagneres  de  Bigorre  in  April  I960. 

(c)  Development  of  Methods  for  determination 
of  the  relative  heights  of  shallow  surface 
features . 

The  aim  of  the  research  carried  out  under  this 


heading  by  Messrs  T.  W.  Rackham  and  E.  D.  Dale,  working 


under  the  supervision  of  Professor  Kopal  at  the  University  of 
Manchester,  has  been  to  study  the  properties  of  lunar  micro¬ 
relief  in  the  region  of  the  maria  by  a  method  based  on  micro¬ 
photometric  tracings  of  the  Pic-du-Midi  lunar  negatives.  This 
micro-relief,  small  as  it  is  when  compared  with  the  topography 
of  the  craters  of  mountain  chains,  attains  some  prominence  on 
the  otherwise  flat  maria.  In  general,,  therefore,  the  study  is 
concerned  with  the  relief  of  the  maria,  and  hopes  to  contribute 
towards  a  clearer  understanding  of  the  nature  of  their  surface. 

Topographic  features  amounting  to  no  more  than  20 
metres  in  height,  and  inclined  at  an  angle  of  the  order  of  one 
degree,  can  be  measured  significantly  by  a  comparison  of  the 
photographic  densities  of  the  negative  on  a  sloping  and  level 
region.  By  a  comparison  of  a  density  of  the  image  of  a  ridge 
(or  dome)  with  equal  density  on  the  mare  at  a  spot  where  the 
Sun  stands  higher  (or  lower)  above  the  horizon,  actual  slopes 
over  the  ridge  (or  dome)  can  be  established.  Graphical  integration 
then  furnishes  a  series  of  spot  heights  over  the  particular  uneveness 
of  the  surface;  and  by  using  a  raster  technique  a  contour  map  of  the 
formation  can  be  drawn. 

Mr.  Dale  applied  this  technique  in  some  detail  to  the 
Serpentine  ridge  in  Mare  Serenitatis,  as  well  as  to  a  representative 
lunar  dome  in  the  Arago  region.  The  quantitative  results  of  his 


work  are  contained  in  an  M.Sc.  thesis,  now  in  preparation, 
to  be  submitted  to  the  University  of  Manchester  in  the  fall 
of  1962;  and  will  also  be  published  as  a  Technical  Scientific 
Note  under  this  contract. 


